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ABSTRACT 

Accelerated  life  testing  of  a  product  under  more  severe 
than  normal  condition  is  commonly  used  to  reduce  test  time  and 
cost.  Data  collected  at  such  accelerated  conditions  is  used  to 
obtain  estimates  of  parameters  of  a  stress  translation  function 
which  is  then  used  to  make  inference  about  the  products  per¬ 
formance  under  normal  conditions.  This  problem  is  considered 
when  the  product  is  a  p  component  series  system  with  Weibull 
*  distributed  component  lifetimes  having  a  common  shape  parameter. 
A  general  stress  translation  function  is  used  and  estimates  of 
model  parameters  are  obtained  for  various  censoring  schemes. 


1. 


INTRODUCTION 


Accelerated  life  testing  of  a  product  is  often  used  to 
obtain  Information  on  its  performance  under  normal  operating 
conditions.  Such  testing  involves  subjecting  test  items  to 
conditions  more  severe  than  encountered  in  the  item's  everyday 
use.  This  results  in  decreasing  the  Item's  mean  life  and  leads 
to  reduced  test  time  and  experimental  costs.  In  engineering 
applications  accelerated  conditions  are  produced  by  testing 
items  at  higher  than  normal  temperature,  load,  voltage,  pressure, 
etc. ,  while  in  biological  applications  accelerated  conditions 
arise  from  large  doses  of  a  chemical  or  radiological  agents.  In 
all  cases  the  data  collected  at  the  high  stresses  is  used  to 
extrapolate  to  some  low  stress  where  testing  is  not  feaslable. 

We  shall  consider  the  problem  of  accelerated  life  testing 
when  the  item  of  interest  is  a  p-component  series  system.  Here 
the  failure  of  any  one  of  p  components  causes  the  system  to  fall. 
An  observation  on  such  a  system  consists  of  the  system  failure 
time  and  knowledge  of  which  component's  failure  caused  the 
system  to  fail.  In  the  case  where  there  is  data  only  at  opera¬ 
ting  conditions,  David  and  Moeschberger  (1978)  describe  some  of 
the  estimation  techniques. 

Several  papers  have  been  written  on  analyzing  accelerated 
life  tests  for  series  systems.  Assuming  that  for  a  given  stress 
V  each  failure  mode  follows  an  independent  log  normal  distribu¬ 
tion  with  parameters  yi(v)  =  V  and  constant  with 

respect  toV,  i  ■  1,  ...,  P,  Nelson  (1973)  obtains  graphical 
estimates  of  and  0^  when  there  is  no  censoring.  Maximum 
likelihood  estimates  of  a^,  8,  and  are  obtained  in  Nelson 
(1974). 

Klein  and  Basu  (1981a)  have  considered  the  above  problem 
when  the  component  lifetimes  are  exponentially  distributed  and 
the  data  is  type  I,  type  II  or  progressively  censored.  Klein 
and  Basu  (1980b)  have  considered  this  problem  when  the  component 


lifetimes  follow  a  Weibull  distribution  with  unequal  shape 
parameters.  In  this  report  we  consider  this  problem  when  the 
components  have  a  common,  unknown  shape  parameter. 

In  section  2  we  present  the  model  we  shall  use  for 
accelerated  lifetests  in  the  competing  risk  framework.  In 
section  3  we  obtain  estimates  of  model  parameters  for  type  II 
and  progressively  censored  data.  In  section  4  an  illustrative 
example  is  presented. 

2 .  THE  MODEL 


The  problem  considered  in  the  sequel  is  as  follows.  Consi¬ 
der  a  p  component  system  with  component  lifetimes  X^,  X2 . 

Xp.  Suppose  that  under  normal  stress  conditions  these  compo¬ 
nents  have  long  lifetimes  making  testing  at  such  conditions 
unfeasible.  To  reduce  test  time  and  cost,  s  stresses,  V^,  ..., 
Vg  are  selected  and  a  life  test  is  conducted  at  constant  appli¬ 
cation  of  the  selected  stress.  We  wish  to  use  this  information 
to  make  inference  about  the  component  lifetimes  under  normal 
stress  conditions. 

Consider  the  following  model  Introduced  by  Klein  and  Basu 
(1980b)  elsewhere. 

At  a  stress  V^,  i  *  1,  ...,  s  assume  that  the  component 
has  a  hazard  rate  given  by 

(x ,  Vi ;  Ctj  ,  ) 


i  *  1, 


gjU,  o^O^,  Bj) 

3  “  l  *  •  •  •  *  p  • 


(2.1) 


For  gj(x,  a)  a  Weibull  form  is  assumed,  that  is 


gj(x,a)  -  at 


a-1 


,  o>  0,  t  >  0, 


(2.2) 


where  a  is  the  same  unknown  constant  for  each  component  and  this 
constant  is  independent  of  the  stress  V. 

For  Aj(V,  Bj)  we  assume  a  model  of  the  form 


(2.3) 


Xij  "  XjtV'  -j*  *  ^  ^  Jo  ^ 


where  0jq(V)  =  1  and  @.,.,(7),  ...,  0<t  (V)  are  k.,  non-decreasing 


Jjl 


jk 


j 


functions  of  V.  The  0j(*)'s  may  differ  from  one  component  to 
another. 

This  model  includes  the  standard  models,  namely,  the  power 

rule  with  (V^  ,  j$  =  3^qV  j  ;  The  Arrhenius  reaction  rate 
model  with  Aj(V,  0^)  =  exp(0jQ  -  0j-^/V);  and  the  Eyring  model 

8,, 

for  a  single  stress  with  Aj(V,  0^)  =  V  J  exp(0j()  -  0j2/V)  as 
special  cases. 


The  model  also  can  be  derived  from  the  interpretation  of 
the  effects  of  a  carcinogen  on  a  cell  as  proposed  by  Armitage 
and  Doll  (1961).  To  produce  cancer  in  a  single  cell,  k  indepen¬ 
dent  events  must  occur.  The  effect  of  an  increased  dose  of  a 
carcinogen  is  to  increase  the  rate  at  which  these  k  events  occur. 
If,  for  the  j1"*1  disease,  this  increase  is  of  the  form 
e*P(3j£0j£(V))  for  1*1,  ...,  kj  the  model  (2.3)  is  obtained. 

If  this  increase  is  assumed  linear  the  model  of  Hartley  and 
Sielkin  (1977)  is  obtained.  Thus  the  model  of  Hartley  and 
Sielkin  is  a  first  order  Taylor  Series  approximation  to  (2.3) 
when  ejA(V)  -  V  for  4  -  1,  ....  kj. 

Consider  an  accelerated  life  test  conducted  at  constant 
applications  of  8  stress  level,  Vj,  ...,  Vg.  Let  X^,  X12, 

•  •  • »  X^  denote  the  component  lifetimes  of  the  p  component 
series  system  put  on  test  at  stress  V^.  Assume  that  the  compo¬ 
nent  lifetimes  are  independent.  We  are  not  allowed  to  observe 
Xu,  ...»  X^  directly  but,  instead,  we  observe 

-  minimum(X^,  ...,  V  and  an  indicator  variable  which 
describes  which  of  the  p  components  is  the  minimum.  We  shall 
use  the  method  of  maximum  likelihood  to  estimate  a  and  ^  m 
(6jo»  •  ••*  3jkj)t  j  *  1 >  ...»  p  for  various  censoring  schemes. 


3. 


ESTIMATION  OF  PARAMETERS 


In  this  section  we  obtain  maximum  likelihood  estimates  for 
a  and  ,  j  “  1,  p  for  type  II  and  progressively  censored 

accelerated  life  tests.  Since  these  estimates  must  be  obtained 
numerically  we  also  show  how  initialestimates  can  be  obtained  by 
a  least  square  technique.  We  then  obtain  point  and  interval 
estimates  for  the  component  parameters  at  the  use  stress. 

3.1  TYPE  II  CENSORING 

For  this  censoring  scheme  n^  systems  are  put  on  test  at 
each  of  the  s  stress  levels  and  testing  continues  until  a 
preassigned  number  r^  have  failed  at  which  time  testing  at  that 
stress  level  is  terminated.  Let  the  ordered  failure  times  of 
the  r^  systems  which  fail  be  denoted  by  Y.^j  ,  Y^)  »  •  •  • » 
j  .  Let  r^  denote  the  number  of  the  r^  systems  which 
failued  due  to  failure  of  the  J1"*1  component  j  *  1,  ...»  p.  The 
contribution  to  the  total  likelihood  from  the  n^^  systems  on  test 
at  stress  is 


p  ri1  ri 

Li  *  nx  exp(-T1(a)Xij)a  1 


‘i 

n  y 

£=1 


i(£) 


a-1 


(3.1.1) 


where 


Ti(a)  -  II  Y®(£)  +  (n1-r1)Y“(r  i  -  1 . s  (3.1.2) 

£*  1  1 


s 

The  overall  likelihood  is  L  *  II  L.  so 

i-1 


-  !  l 

i-1  J-l  3 


( 


So 


-  Tt(a)exp 


So  Vjt<V 


+  r.  +  (a-l)P 


(3.1.3) 


where 


■ 


P1- 

£=1 

The  likelihood  equations  are 

HZrm  X  Vj*(V»  +  ir+  pi  - 0 


i=l 


j-1  n=i 


(3.1. A) 


where 


i1)<a)  '  *?<«  ™1W  +  <ni  -  ri,r?  **!<>) 


i  “  1,  . . . ,  s,  and 


^  -  °  -  I  9jg<*1)i»lj-t1(a).xP(  r  Vj»(V» 

JU  i—i  jC—[) 


j  a  1»  • • • 9  P  u  If  . . . ,  k 


j 


(3.1.5) 


The  solutions,  a,  Bj^,  j  -  1,  . ..,  p,  l  «  1,  ....  kj  of  3.1.4 
and  3.1.5  are  the  respective  maximum  likelihood  estimators. 
The  second  partial  derivative  of  L  are 


62  •**!. 


S  IT 

£  A.T^(a)  +  l  —where  (3.1.6) 

_  v  1  1  j  i  _.2 


6a  i=l 


i=l  a 


and, 


o  9nL  ,  , 

-  Te~ 68  =  0  i  *  8  i  “  1.  •••»  p  (3.1.9) 

ju  Wgw 

8-1 . .  "■° . kj  ”-° . y 

To  find  E(Tita)) ,  E(t£1)(ci)),  E(T<2)(cO)  note  thet  the 
unordered  Y^'s  have  a  Weibull  distribution  with  density 

f(yiA)  -  «VS5le*P(-V?i>;  yU  >  °*  1  =  1 . .  • 

The  mth  order  statistic  from  a  Weibull  has  density 

£(W  ‘  rnfLi  ‘“VO1-  yi(m)]"-1 


expf-^-w-DX^.  .),  0  <  y  <  «. 


(3.1.10) 


E<1?(«)>  ■  /  TS-TKVO!  >aK>1°)ya-1m-«xp(-y“x1))'°-1 

exp(-(ni-nH-l)y0tA1)dy 
n±!  00  u  m  1 

=  05ri)!(n1Mn)!Xi  /  ul1-e"u]m"  exp(-(ni-^l)u)du 

(n.lm-1  .  fm-l)  00 

.  -  J  Jo  (-l)  U  J  l 

_  n.  m-1  .  [m-1  , 

-  -3_  i  Y  1 

xi  i"Jk-o  k  * 

.  4 


Hence 


i  -nrt-l+k)  * 


i 


Let  a  -  -E( 


~~*'7i — ) ,  Let  C.  be  the  column  vector 
6a2  j 


V 


jo 


,62^L  . 


>  j  =  •••»  P»  snd 


let  Bj  be  the  k^+1  by  k^+l  matrix  whose  elements  are 

.  62  m  L  x  -  , 

bu+l,w+l  e^'<5$.  6ft  *  u  “  °»  •**»  kj»  w  =  °»  ••• 

ju  jw  J 

The  inverse  of  the  asymptotic  correlation  matrix  is 


‘j  ‘ 


,-l 


B  C 

B^  0 

CT  a 

where  B  = 

B2  . 

0  "B 

P 

.  . 

(3.1.14) 


and  CT  =  (cT,  C~,  ...,  C^.  The  value  of  £  ^  is  obtained  by 
substituting  3.1.11,  3.1.12,  and  3.1.13  in  the  appropriate  places 
in  3.1.6,  3.1.7,  and  3.1.8. 

By  a  theorem  of  Rao  (1973) 


l  = 


-1  -1  T 

B  +FE  ■LFi 


-1  T 
-E  T 


-FE 

„-l 


-1 


where 


F  -  B_1C,  and  E  =  a-CTB_1C. 


(3.1.15) 


A  consistent  estimator  of  I  is  obtained  by  using  a  and  8^ , 

j  ■  1 .  p  in  the  appropriate  expressions  3.1.11,  3.1.12, 

and  3.1.13. 


3.2  TYPE  I  PROGRESSIVE  CENSORING 


For  this  censoring  scheme  N 
stress  Let  T^,  •  •  • ,  T 
At  censoring  time  x^,  l  -  1,  .. 
items  are  removed  from  the  study 


^  items  are  put  on  test  at 

. u  be  fixed  censoring  times. 
1Mi 

•>M^-1  ,  a  fixed  number  c^, 

.  At  time  x.u  either  a  fixed 


number  items  are  removed  from  the  test  and  testing 

continues  until  all  remaining  items  have  failed  or  the  test  is 

terminated  with  a  random  number  c^M  items  still  functioning. 

i 

This  test  scheme  has  the  advantage  of  allowing  for  some  items 
with  extremely  long  lifetimes  to  be  encluded  in  the  study. 
Clearly  the  usual  type  1  fixed  time  censoring  is  a  special  case 
of  this  consoring  scheme  with  =  i  and  c.^  random.  We  assume 
that  N±  is  sufficiently  large  so  that  at  least  c^,  £  =  1,  . .., 
M.1,  items  are  still  functioning  to  be  censored. 

Mi 

Let  n  =  N.  -  £  c..  be  the  number  of  systems  which 

i=l  1 


failed  and  let  Y^,  Y^.  •••,  Yin  denote  the  failure  times. 

Let  r^  denote  the  number  of  systems  which  fail  from  cause  j  at 
stress  V^,  j  =  1,  ....  p,  i  =  1,  . ..,  s.  The  items  on  test 
at  stress  i  =  1,  ....  s  contribute 


L 


i 


ss 


p 

n 


rn  ni 

exp^X^T^Ct^ot  1 


n 

£=1 


-1 


i£ 


(3.2.1) 


to  the  total  likelihood  where 


Ti(a) 


l  * 


£=1 


a 

i£ 


M. 


£**1 


a 

Ti£ci£ 


(3.2.2) 


Again  the  total  log  likelihood  is  given  by  3.1.3  with  T^a)  as 


in  3.2.2  and  P  ■  ][  tn  Y.  .  Maximum  likelihood  estimators  a 

1  £-1  1 

A 

and  ,  j  ■»  1,  ....  p  of  a  and  are  obtained  by  solving 
numerically  3.1. A  and  3.1.5  with 


i  ni 

Ti  ^  (a)  “  Yi£  Mi£  +  Ti£  37rTi£ci£* 


(3.2.3) 


The  second  partial  derivatives  of  L  are  given  by  3.1.6,  3.1.7, 


3.1.8,  and  3.1.9  with 


ni  Mi 

Ti  ^  ^  Yi£^Yi£^  +  cliTll('2nTil) 


(3.2.4) 


We  now  calculate  E(T^(a)),  E(T^(a)),  and  E(T^(a)). 

Consider  any  of  the  s-stress  levels.  For  notational  con¬ 
venience  we  shall  suppress  the  subscript  i.  Let  N  items  be  put 
on  test.  Let  0  =  tq  <  ...  <  tm  <  TMfl  =  °o  be  censoring  points. 

£,  1  =  1,  ...,  M,  c^  items  are  removed  from  test.  For 


At  time  Tn>  £ 


£  =  1,  ....  M-l,  c^  is  a  fixed  constant  while  c^  is  either  ran¬ 
dom  or  fixed  depending  on  if  testing  is  terminated  at  or  not. 


Let  n 


fail. 


M 

N-  l 

£=1 


c^  be  the  number  of  items  which  are  observed  to 


Let  Y^,  ....  denote  the  failure  times  of  the  n  failures. 
For  £  ■  1,  ...,  n,  Y^  =  min(X^£,  ...,  V  where  X^£  is  the 
failure  time  of  the  jbh  component  of  the  £th  item  which  fails. 

By  assumption  the  Y^'s  have  a  Weibull  distribution  with  survival 
function 


F(y)  *  exp(-Xya),  y  >  0,  a,  X  >  0  (3.2.5) 

Let  f^  denote  the  number  of  failures  in  the  interval 

[T£_lf  Tg) ,  i  *  1,  ....  M+l.  Let  U£k>  k  "  1 . f£’ 

£  *  1 . M+l  denote  the  failure  times  of  those  f^  items 

which  fail  in  the  interval.  Let 


F£  -  P (Y  <  Tjl)  -  1  -  exp(-XxJ)  i  “  1 . M+l  (3.2.6) 

denote  the  probability  an.  item  fails  before  time  T^,  and  let 

'  1  '  Fr 

Cohen  (1963)  shows  that  for  a  fixed  constant 


E(f*) 


NF, 


for  £ 


£-1  c. 


(N  -  l  =£)  (Ff  -  Ff  .)  for  £  -  1,  ....  M+l 
k-1  Fv 


(3.2.7) 


For  random  3,2,7  holds  for  A  ■  1,  . M, 

fA 

Now  for  A  -  1,  M  let  S  -  Y  U?,  and  note  that 

j-1 

f  A  ,,a  i,  , 

E(S  )  =  EE(  l  UAj|1:A; 

j-1 

-  "V'KjI'm  <  »tl  <  V  • 

The  conditional  density  of  given  t£)  is 


W"> 


{  aXuQ  ^exp(-Xua) 

fa-fa-i 


lf  tA-1  <  u  <  ta 


otherwise 


Thus 


80 


A-l  c 


E(SA)  "  X(N  "  ^  Ii)((KTA-l+1)FA-l  '  (xtJ+1)F£)A  =2 . M 


j=l  F 


(3.2.8) 


and 


E(S£)  -  j  N(1  -  (XtJ+l)F1) ,  A-l 


We  now  consider  separately  the  two  possible  cases  for  c^. 
Case  1,  fixed. 

For  this  case  note  that  3. A. 8  holds  for  A  -  M  +  1.  So 
M+l  M+l  _  ' 

xJi  E<st>  -  » 


M  c£  M+l 


-  I  “  I  KXtJ,+1)F.  ,  -  (Xr.+1)F. )  . 
A-l  F.  J-1  2  J  2  2 


since  tq  ■  0,  Fq  -  1  and  F, 


M+l 


0  we  have 


M+l  .  M 

I  E(St)  ■  $  (»  -  I  CjCXtJh))  . 
&=1  A=1 


»i  -  l 


Thus  E(Ti(a))  = 


t  i  *  1 1  i •  •  f  Si 


(3.2.9) 


Similarly, 


EClJ^Co))  =  {N±(l  -  Y  -  07!  Xp 


J.  w 

"  3.  Ci£[,  /a  u  3711  exP<"u  “  XiTi£)du 

X,=l  a1tiA 

-  finX^  -  X^T^ln  X^T^]  ,  where  y  Is  Euler's  constant. 

The  integral  must  be  evaluated  numerically. 

Also  by  similar  computations 

E(T<2)(ci))  -  Nt(^-  -  2Y  +  y2)  -  22nX1(l-y)  +  (^X±)2 


i  -  -(u~XlTiO> 

-  I  cij.l  /  (2nu)2ue  1X1  du 

‘-1  Vh 


-  2  2m\ 


00 

{y 

AiTi  l 


(  0nu)ue 


-(u-X.t,„) 


du  +  X. 


XiT±£<0n  V«>  (  (Xi/TU>  >  1 } 


Case  II. 


random. 


For  this  scheme  all  testing  stops  at  time  with  cM  items 
still  functioning.  Clearly 


So  from  3. A. 7 


M-l  M 

Cm  “  n  ~  I  c £  I 

£-1  £=1  * 


r  c£ 

,{*  ’  F«<N  -  Ji  i? 


(3.2.12) 


E  Is*  ■  ^fi*a-(»TS+i)pM  -  V  ^  katJ+df,  -  (at“+dfmi) 

.:L=1  ^“1  E„ 


Hence 


M±_1  F 

E(T  (a))  -  -f^N  ±F  -  l  ryd-^)}  (3.2.13) 

i  1  £=1  1Jt  F 

“i 


Similarly, 


,n  .  i  1M 

E<Tl(“»1^Nlf/  ‘ufeiue^du-  M^+f^t^  *0^)1 

M±-l  (  ViM. 


-  AiT-j  WAjtJp+AjT^  MA^^F^/Fjj  ),  (3.4.14) 


.Ot  xzr  /- 


X  Ta 
iTiM 


/o\  -1  «**  X  iXl. 

E(T}  }(a))—±-  N  [/  u(  0nu) 2e_udu-2  3»X./  u  3nue"udu 

■*■  1  1  n  a  . 


\a‘  *  0 


+<*  «"  Xl>  *«/> 


"i"1  Vi*  .fc  ,1 


11  Va 


e  udu-2®n  X .  /  ue"udu 

i  >  _a 


Vu 


(3.2.15) 


+<fa  A±)  2 U-F^/F^+A^KM^  2-< 


ra  >21 

TU'  J 


The  asymptotic  covariance  matrix  is  obtained  by  making  the 
appropriate  substitution  in  3.1.6,  3.1.7  and  3.1.8. 


3.3  INITIAL  SOLUTIONS  TO  LIKELIHOOD  EQUATIONS 

To  solve  the  likelihood  equation  numerically  initial 
estimates  of  the  parameters  which  are  close  to  the  maximum 
likelihood  estimates  are  needed.  To  obtain  such  estimates  we 
shall  first  obtain  an  estimate  of  a,  then  transform  the  data 
to  exponential  observations  and  apply  a  least  squares  technique. 

Consider  any  one  of  the  s  stress  levels,  v^.  Let  be  an 
estimate  of  a  based  on  observations  at  this  stress  level  only. 
Such  estimates  can  be  obtained  by  using  techniques  described  in 
Mann,  Schaffer  and  Singpurwalla  (1974)  or  by  graphical  methods 
described  in  Nelson  (1972) .  These  estimates  are  then  pooled  to 
obtain  an  estimate  a  of  a.  If  the  a^s  differ  too  much  from 
one  stress  level  to  the  next  this  will  cast  doubt  on  the 
assumption  of  equal  shape  parameters. 

To  obtain  estimators  of  the  j5's  we  first  make  the  trans¬ 
formation  =  Y^a,  i*  1,  ...,  s,  i,  =  1,  ...»  n^.  If  a  is 
equal  to  the  true  a  then  will  have  an  exponential  distri¬ 
bution  with  hazard  rate  A , . 

%  1 

Let  T^(ct)  be  defined  by  3.1.1  or  3.2.1.  One  can  show 
that,  using  information  collected  at  a  single  stress  level  only, 
the  maximum  likelihood  estimator  of  A.,  is 

^rij 


say.  For  sufficiently  large  n 


E<y  *  Wy  -  Jo 


(3.3.1) 


(3.3.2) 


and  Var(n^)  -  ^1;jE(T1(o)) 


Least  squares  estimates  of  the  0^*8  are  obtained  by  using 
weighted  least  squares  as  described  in  Draper  and  Smith  (1966, 

A 

p.  77-81).  Let  ri  be  the  sXl  column  vector  of  the  r^’s.  Let  B 
be  the  s  x  +  1  matrix  defined  by 

6  =  (0jA<Vi>>»  1  m  0 . kj  1=1 . 8 


and  let  B  be  the  column  vector  (Bjq» 
interest  is 


3.,  )T.  The  model  of 

1 i 


n  =  SB  +  e  (3.3.3) 

where  E(e)  ■  0  and  Var(e)  ■  A  where  A  is  the  diagonal  matrix 
with  elements  Var(ri^)  along  the  main  diagonal.  The  weighted 
least  squares  estimators  of  are 

ij  -  (0ta* 1 6) _1 eTA-1n  •  (3.3.4) 

The  variance  of  these  estimators  is 

Var ( 3^ )  -  (0TA-10)~l  .  (3.3.5) 


3.4.  ESTIMATION  OF  USE  STRESS  PARAMETER 

Suppose  an  accelerated  lifetest  as  been  conducted  according 

to  one  of  the  censoring  schemes  discussed  in  sections  3.1  or 

3.2.  Let  jij  ■  (Bq . 8^  )  and  a  he  the  maximum  likelihood 

J  A 

estimators  of  and  a  respectively.  Let  ^jj(^jj)  he  the 
(estimated)  asympotic  variance-covariance  matrix  of  .  Let 
Eja(rja)  be  the  kj  +  1  *  1  column  vector  whose  8.  element  is  the 
(estimated)  asymptotic  covariance  of  0..  and  a,  £.  ■  0,  ...,  k. , 

.  ^  j*  j 

j  -  1,  ...»  p.  Let  ^^((O^jj)  he  the  (estimated)  asymptotic 
variance  of  a.  Let  0  be  the  matrix  whose  elements  are  all  zero, 

A  A  A 

then  the  asymptotic  variance  matrix  of  (j^,  B^,  •••»  a)  is 
of  the  form 


(3.4.1) 


Z  Z 
pE  « 

Z  T  a  2 
pa  aa 


Let  Z  denote  the  estimated  value  of  Z. 


We  shall  use  this  Information  to  make  Inference  about  component 

life  under  some  use  stress,  V  . 

u 

Recall  that  the  scale  parameter  of  the  time  to  failure 
distribution  Is  given  by 

\>  1 

Aju  ■  exp  l  £0^ £<Vu)  ,  j  *  1 . .  .  (3.4.2) 


at  the  use  stress.  The  maximum  likelihood  estimator  of  Xju  is 
given  by 


“  exp  £  8j£®jji(Vu)  ,  j  —  1,  ...,p.  (3.4.. 


This  estimator  is  biased  by  a  factor  of  exp(o2u/2)  where 


"Jk  vru'“jjvx’  WJ1'V 


jkj  '  u''  * 
(3.4.4) 


;An  unbiased  estimator  of  is  given  by 

Xju  “  *ju  exP(“°ju^2^»  J  “  1 . .  •  (3.4.5) 

Asymptotic  (1  -  y)  x  100*  confidence  intervals  for  Aju  are  given 

^juexp(”Zl-Y/2°ju)’  ^juexp(Zl-y/2aju)  (3.4.6) 

where  zi-y/2  t*'e  l-Y/***1  percentage  point  of  a  standard 
normal  random  variable. 


Consider  now  the  cumraulative  hazard  rate  of  the  time  to 
failure  distribution  of  the  j**1  component  which  is,  at  time 
t  >  0  and  stress  Vu,  given  by 


a 


AJu(t)  =  t  exp 


Jo  WV 


j  -  1 . .  (3.4.7) 


The  maximum  likelihood  estimator  of  Aju(t)  is,  by  the  invari¬ 
ance  property. 


AjU(t)  “  t°exp 


Now,  asymptotically 


i  vv<v 


Aju(t)  and  variance 


Ju 


,  j  —  1,  ...»  P*  (3.4.8) 

(t)  has  a  normal  distribution  with  mean 


oJ„(t)  -  (1,  eji(vu) . e11r  (Vit),  *.(t)) 


JU 


jkj  u 


J  "  1.  •••»  P,  t  >  0 


zjj  Zja 

E  T  E  2 
ja  aot 


'9jl  <Vl 


%<V 

<bn  (t) 


(3.4.9) 

A 

Hence  Aju(t)  is  a  biased  estimator  of  A^u(t).  A  reduced  bias 
estimator  of  A^u(t)  is  given  by 

*\j  %  A 

Aju(t)  -  AJu(t)exp(aju(t)/2)  j  -  1 . P,  t  >  0  (3.4.10) 

which  also  has  reduced  mean  squared  error 
A  (l-y)x100Z  confidence  interval  for  Aju(t)  is  given  by 

(^ju(t)exp(‘Zl-Y/2°ju(t)) *  ^ju(t)exP(Zl-y/2aju(t)))  * 

Let  F  (t)  “  exp(-Aju(t)) ,  t  >  0  be  the  survival  function 
of  the  jth  component.  The  maximum  likelihood  estimator  of 

Vt}  18 

Fju(t)  "  exp(-AJu(t))  J  -  1, 


•  e  • 


p,  t  >  0 


(3.4.12) 


Approximate  (1~y)x100%  confidence  intervals  for  F^u(t)  are  given 
by 


<V° 


e*P<zl-r/2°Ju>  f  ,°*H-Zl-y/2°)u\ 

,  Fju(t)  ), 


1 »  ••  • I  p 

(3.4.13) 


3.5. 


DEPENDENT  RISKS 


The  assumption  of  Independent  causes  of  failure  may  be 
relaxed  to  include  a  class  of  Weibull  distributions  with  depen¬ 
dent  causes  of  failure.  To  illustrate  how  this  may  be  done  we 
shall  consider  the  bivariate  case  with  the  obvious  extension  to 
more  than  two  risks. 

Let  X^,  X2  denote  the  failure  times  of  the  two  components 
in  a  series  system.  Let  be  the  time  until  the  system  fails 
due  to  failure  of  the  first  component  alone,  U2  the  time  until 
failure  from  the  second  component  alone,  and,  U^2  the  time 
until  simultaneous  failure  of  both  components.  At  a  stress  V 
assume  that  U^,  U2,  Uj,2  are  independent  Weibull  random  variables 
with  shape  parameter  a  and  scale  parameter  Aj(V,  3j)»  j  ■  1,  2, 
12,  given  by  2.1.2.  Clearly  X^  =  min(U^,  U^2)  and  X2  •  min 
(u2,  u12).  X^  and  X2  are  both  distributed  Weibull  with  shape 
parameter  a  and  scale  parameter  ^(V,  B^)  +  X12(V,  012)  , 

X2(V,  $2)  +  X12(v,  —12^’  respectively.  The  joint  survival 
function  of  (X^  X2)  is  given  by 

F(xlt  x2)  -  expC-Aj^V,  B1)x1-X2(V,  612)max(x1,x2)a) . 


This  is  the  bivariate  Weibull  distribution  proposed  by  Lee  and 
Thompson  (1974) . 

To  estimate  the  parameters  we  perform  an  accelerated  life 
test  as  described  in  the  previous  sections.  The  failure  causes 
are  now  failure  from  the  first  component  alone,  the  second 
component  alone,  and,  simultaneous  failure  from  both  components. 
Estimators  of  jj 2 ,  0^2 ,  and  a  can  be  obtained  as  before. 

Estimates  of  the  parameters  of  X^  and  X2  now  follow  in  a 
straight  forward  manner. 


EXAMPLE 


As  an  example  of  these  procedures  we  shall  consider  an 
example  given  In  Nelson  (1974a) .  The  problem  Is  to  analyze  an 
accelerated  life  test  conducted  on  Class-H  insulation  systems 
for  electric  motors.  There  are  three  possible  types  of  insula¬ 
tion  failures  corresponding  to  distinct  parts  of  the  insulation 
system,  namely  Turn,  Phase,  and  Ground.  The  failure  cause  is 
determined  by  an  engineering  examination  of  the  failed  motor. 

The  purpose  of  the  experiment  is  to  estimate  the  average 
life  of  such  insulation  systems  at  a  design  temperature  of 
180°C.  A  median  life  of  20,000  hours  is  necessary  for  the 
satisfactory  performance  of  these  insulation  systems.  To 
reduce  test  time  and  cost  an  accelerated  life  test  was  conducted 
at  4  accelerated  temperatures,  namely,  190°  C,  220°  C,  240°  C, 
and  260°  C. 

The  accelerated  life  test  was  conducted  by  putting  10 
motors  on  test  at  each  of  the  4  stress  levels.  Motors  were  run 
until  they  failed,  then  the  cause  of  failure  was  found  and 
Isolated  and  motors  were  run  until  a  second  failure  occurred. 

The  results  of  this  study  are  reported  in  Nelson  (1974a) .  The 
data  followed  a  log^Q  normal  distribution  so  the  Weibull  theory 
results  do  not  apply. 

To  illustrate  the  results  of  the  previous  section  Nelson's 
example  is  reporduced  by  simulating  the  life  test  using  a 
Weibull  model  with  shape  parameter  1  for  each  failure  cause. 

The  shift  parameters  are  chosen  by  fitting  an  Arrhenius 
Reaction  Rate  model  to  the  estimated  component  medians  obtained 
by  Nelson.  The  model  is 

Xt(V;  gj)  -  exp(Bj0  +  Bj^iOO),  j  -  1,  2,  3  (4.1) 

where  0^(V)  ■  -1000/V  for  J  ■  1,  2,  3  and  V  is  the  temperature 
in  degrees  absolute.  The  absolute  temperature  is  273.16  plus 
the  centrigrade  temperature.  The  values  of  (Bjq*  6j.j),  j  ■  1, 


2,  3  are  as  follows; 

Table  4.2  True  Values  of  Bq,  B^ 

g0 _ gl 

Turn  8,2607  8.0106 

Phase  3.7748  6.1253 

Ground  13.0340  10.6487 

Twenty  Weibull  observations  were  generated  at  each  of  the  four 
stress  levels.  The  data  are  In  Table  4.1. 

Using  the  data  at  each  stress  level  only  the  two  order 
statistic  estimates  of  ot  discussed  in  Derbey  (1966)  are 

AAA  A 

=  2.224,  0.2  =  1.103,  =  1.260,  and  ot^  =  1.155.  This 

*\j 

suggests  an  initial  estimate  of  a  *  1.425  for  a. 

The  data  is  now  transformed  by  letting 
i  =  1,  ...,  s,  SL  =  1,  ...,  n^  and  the  least  squares  procedure 
of  section  3.3  is  applied.  The  resulting  initial  estimate  of 
the  B's  are  as  follows 

Bturjj  *  (4.3136,  4.9144),  BpHASE  -  (1.5206,  3.8801)  , 
and 

^GROUND  "  <6'1364’  6-1194)  * 

The  maximum  likelihood  estimates  are  obtained  by  using  a 
two  stage  Newton-Raphson  procedure.  In  the  first  stage  the 
likelihood  in  maximized  with  respect  to  a  using  the  B's 
obtained  in  the  previous  stage.  In  the  second  stage,  using 
this  a,  the  likelihood  is  maximized  with  respect  to  the  B's. 
The  procedure  is  terminated  when  the  relative  increase  in  the 
likelihood  is  no  more  than  .0001. 

In  this  case  the  two  stage  procedure  terminated  after  19 

A 

steps.  The  maximum  likelihood  estimators  are  a  -  1.0995, 

Irma,  -  (6.6142,  7.5033),  Iphase  -  (4.1230,  6.6608),  and 
^GROUND  *  (8*2946,  8.6483).  The  estimated  covariance  matrix 

of  (£rURN»  IpHASE’  ^GROUND*  Q) ’  co»Puted  fr°“  the  results  of 
section  3.1  with  r^  ■  n^,  is 
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At  the  use  stress  of  180°C  the  estimates  of  component 
survival  at  a  mission  time  of  20,000  hours  are  .0763  for 
turn  failures.  90%  confidence  intervals  for  components 
survival  at  20,000  hours  are 

(.0013,  .3688)  for  turn  failures, 

(.0004,  .7872)  for  phase  failures,  and 
(.00242,  .8168)  for  ground  failures. 

The  maximum  likelihood  estimates  of  the  scale  parameters 
at  180°C  are 

Aturn  *  .0000480,  *pHASE  “  .0000255,  and  *GR0UND  -  .0000206. 

The  reduced  bias  estimates  are 

*TURN  "  'O000*54*  ^pHASE  “  .0000221,  and  *GR0UND  =  *0000184. 

90%  confidence  intervals  for  the  shape  parameters  are 
(.0000275,  .0000839)  for  turn  failures, 

(.0000105,  .0000619)  for  phase  failures,  and 
(.0000093,  .0000455)  for  ground  failures. 
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